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FIG. 1. The diagram of the two-photon process e + e —> e + e + R. 
1. INTRODUCTION 

The purpose of this work is to develop an efficient event generator for the process of the 
two-photon resonance production e + e~ — > e + e~R in the so-called single-tag mode, when one 
of the final electrons^ is scattered at a large angle and detected. Such generator is needed 
for simulation of experiments on the measurement of the meson-photon transition form 
factors. The generator GGRESRC described in this work was used for the measurement of 
the transition form factors for the ir°, rj, r]', and r\ c mesons with the BABAR detector. To 
achieve required accuracy (~ 1%), the radiative corrections to the Born cross section are 
taken into account. In particular, extra photon emission from the initial and final states are 
simulated. 

In the two-photon process e + e~ —> e + e~R, the virtual photons, radiated by the colliding 
electrons, form a C-even resonance with the four- momentum k = k\ + ki (see Fig. [1]). Let 
Q\ be the absolute value of the four-momentum squared, carried by the space-like photon 
connected with the tagged (detected) electron , while Q\ be the same parameter for the 
untagged (undetected) electron (Q\ 0). The transition form factor is determined from 
the measured differential cross section (dcr/dQDdata and the MC calculated cross section 
(dcr/dQ^)Mc: 

I cidata //o2\|2 (da/dQ^data i r MC ( rfi\\1 f1\ 

IW^)! " (da/dQlW 1 F ^ r{Q2)1 ' (1) 
where \FM^ R (Ql)\ 2 is the transition form factor used in MC simulation. 



Unless otherwise specified, we use the term "electron" for cither an electron or a positron. 
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2. BORN CROSS SECTION 

To describe the process e + e~ — > e + e~R we use the notations defined in Fig. [TJ and the 
following six invariants: 

t\ = Q 2 = k 2 } t 2 = Q\ = ^2) 

s 1 = (p[ + k) 2 , S2 = (p' 2 + k)\ (2) 
s = (Pi+P2) 2 , W 2 = k 2 = (k l + k 2 ) 2 . 

The differential cross-section for this process in the lowest QED order is given by [l|: 



a 2 / (W^-frfr dV,dVg 
a 167r%t 2 \/ (pip 2 ) 2 -mf ^ £ 2 ' 



where a is the fine structure constant, m e is the electron mass, E[ (i=l,2) are the energies 
of the scattered electrons and 

2j — 4p x p 2 (Ttt + ^Pi P 2 <^T5 + -^Pi P 2 a ST + Pi P 2 OSS (4 J 

+ 2|p^ _ p^ _ |rTT cos 20 — 8 Ip^p^ 1 tt5 cos 0. 

Here is the angle between the electron and positron scattering planes in the center-of-mass 
(cm.) frame of the virtual photons, a a b are the 7*7* — > R cross sections for unpolarized 
transverse (a, b = T) and scalar (a, b = S) photons. The interference terms containing the 
functions t^ arise due to virtual photon polarization. The function ttt is the difference 
between cross sections for transverse photons with the parallel and orthogonal linear polar- 
izations: Ttt — en — while the cross section for unpolarized photons is ott = +cr±)/2. 

The effects of the strong interaction are completely contained in the functions a a ^ and 
T ao . All other functions entering in Eq. (jlj) are calculable with QED. The expressions for 
the virtual photon density matrices pf + , pf~, pf°, p°° (i = 1,2) can be found in Ref. 

In the case of the pseudoscalar meson production, only the functions ott and ttt are 
non-zero, and 

Ttt — — 2,&tt a. 

The cross section &tt for a narrow pseudoscalar meson 
with the mass Mr can be written in term of the transition form factor: 



ivi , 



R 



F(QIQ 2 2 ) 



F(0,0) 

where T 77 is the meson two-photon width. It should be noted that some two-photon event 
generators neglect the term with ttt- This approach may be reasonable only for study of 
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two-photon processes in the no-tag mode, when both the electrons are scattered at small 
angles. The t?t term gives a sizable contribution to the differential cross section da /dQ\ at 
large Q\ and should be taken into account in simulation of single-tag experiments. 

In the GGRESRC events generator we perform integration of the differential cross section 
using invariant variables fl2]). For a narrow pseudoscalar resonance, Eq. ([3]) can be rewritten: 



4a 2 r 77 F{t x ,t 2 

da — 



F(0,0) 



di 2 dtidsids 2 , . 



where A^s, S\, s 2 ,ti,t 2 , M^rnf) is the Gram determinant Q]. The physical region in the 
variables si, s 2 , *i, t 2 is defined by the condition A 4 < 0. The function B coincides, up to a 
factor, with the function S (Eq. fll])) for pseudoscalar mesons. It was calculated in Ref. {4] 
and is given by 



B = U 1 t 2 B 1 -AB 2 2 +m 2 e B^ (7) 



where 



"c 1 



B x = (A Pl p 2 - 2 Pl k 2 - 2p 2 k x + k x k 2 ) 2 + (hk 2 ) 2 - 16tit 2 - 16m* 

b 2 = ( P ip 2 )(hk 2 ) - (pih^h), (8) 

B 3 = t x (2p x k 2 - k x k 2 f + t 2 {2p 2 k x - hk 2 f + AmKhhf, 

To describe the Q\ and Q\ dependencies of the transition form factor F(Qf, Q 2 ), the two 
options are implemented in the generator: F(Q 2 ,Q 2 ) = F(0, 0), and the vector-dominance 
model (VDM) 

1 (9) 



' ' (l+Qf/A2)2(l + g2 /A2) 2' 

where A = m p for ir°, T], 7/ production, A = mj/^ for T] c , and A = mx for 775. The Q 2 
dependence of the \F\ 2 calculated with Eq. ()9]) at Q\ = for A = m p , is shown in Fig. |2j 

Four-dimensional Monte-Carlo integration of Eq. ([6|) is performed using the method de- 
veloped for the GALUGA two-photon event generator |5( . In this method, in particular, the 
invariant variables are generated in the order t 2 , t\, s\, s 2 . This allows to set a restriction 
on Q\ at the beginning of the event generation and significantly increase the generation 
efficiency for single-tag events. The values of the generated invariants, are then used to- 
gether with a random azimuthal angle of the system of the final particles to calculate the 
4-momenta of the scattered electron, positron, and produced resonance. The formulae to 
do this can be found in Ref. J(5|. The main decay modes for 7r°, 77, and 77' are also simulated 
according to Ref. 0. 
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FIG. 2. The Q 2 dependence of the form factor |F| 2 at Q\ = 0, A = m p =0.7755 GeV. 

The total widths of the rj c and T/& resonances are comparable or even larger than the mass 
resolution of modern detectors. Therefore, the mass distributions for these resonances are 
generated using Breit-Wigner distributions. 



3. RADIATIVE CORRECTION 



In the no-tag mode, when both the electron and the positron are scattered predominantly 
at small angles, the radiative correction to the Born cross section is expected to be small, 
less than 1% [7|. The situation changes drastically in the single-tag mode, at a large electron 
scattering angle. At large Q 2 the correction due to extra photon emission from the initial 
state may reach several percents and should be taken into account in simulation. 

The process-independent formula for the radiative correction in the next-to-leading order 
for two-photon processes in the single-tag mode was obtained in Ref. js]. The main contri- 
bution to the correction comes from the vertex of the tagged electron. The corresponding 
contribution of the untagged-electron vertex is expected to be smaller than 0.5% and ne- 
glected. Fig. [3] shows the diagrams taken into account in Ref. |8|]. They substitutes for the 
left-hand vertex in Fig. [TJ 

The cross section for a single-tag experiment is given by: 

da = da B (l + 5) = da B (l + 5' + 5 VP ), (10) 
where da B is the lowest-order cross section for the two-photon process given, for example, 
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FIG. 3. Diagrams used for calculation of the radiative correction. 

by Eq. The total radiative correction is separated into two parts: 

i. 5', which includes the virtual correction due to the interference between the diagrams 
(a) and (c), soft-photon part of diagrams (d) + (e), and the corrections due to real 
photon emission from the initial (diagram (e)) and final (diagram (d)) states, 

ii. 5yp, the vacuum polarization correction due to the interference between the diagrams 
(a) and (b). 



To obtain 5' we have used the result of Ref. [8( for the total radiative correction, removing 



from it the contribution of the vacuum polarization diagram, 5 e (in Ref. jsj only electron 
contribution was taken into account). The resulting 5' is given by 



In 1 3 



+ (11) 



^max 4 

where r max (<C 1) is the maximum energy of the photon emitted from the initial state in 
units of the beam energy E^, L = In (Q 2 /m 2 ), and Q 2 is the absolute value of the momentum 
transfer squared to the electron. The formula does not contain any restriction on the energy 
of the photon emitted from the final state, i.e. the cross section given by Eq. ([TO]) is calculated 
for the case when the tagged electron is allowed to radiate a photon of any possible energy. 
The values of the correction 5' for nine representative sets of Q 2 and r max are listed in 
Table |TJ In the Q 2 region from 1 to 100 GeV 2 available for experiments at 5-factories, the 
correction reach 5-7% even with the relatively loose restriction (r max = 0.1) on the scaled 
energy of the undetected photon emitted from the initial state. 

The correction 5' is partly compensated by the vacuum polarization correction Syp, for 



which we use the results of Ref. 



10j, which includes the contributions from the e, ji, r 
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TABLE I. The correction 5' (%) for the various values of r max and Q 2 . 



Q 2 (GeV 2 ) 


r max — 0.03 


r max — 0.05 


fmax — 0-1 


1 


-9.1 


-7.4 


-5.2 


10 


-10.6 


-8.6 


-6.0 


100 


-12.1 


-9.8 


-6.8 



c 

o 



o 
o 



— Radiative correction due to vacuum polarization |- 
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FIG. 4. The vacuum polarization correction as a function of Q 2 . The curve "All" shows 5yp calcu- 
lated in Ref. [10] with account of contributions from e, /i, r, and hadrons. The curve "Electrons" 
- represents the contribution only from electrons, 5 e . 

leptons, and hadrons. The Q 2 dependence of Syp is shown in Fig. H] in comparison with S e . 
The values of the total correction 5 = 5' + 5yp calculated for for nine representative sets of 
Q 2 and r max are listed in Table [TT1 



TABLE II. Total radiative correction 5 = 5' + 5yp- 



Q 2 (GeV 2 ) 


r max — 0.03 


fmax — 0.05 


fmax — 0.1 


1 


-5.9 


-4.3 


-2.0 


10 


-5.6 


-3.7 


-1.0 


100 


-4.8 


-2.6 


+0.4 



The emission of the hard photon by the electron distorts the kinematics of two-photon 
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event. To model how this effect influences the detection efficiency, the event generator 
includes generation of extra photons emitted from the initial and final states. 

3.1. Simulation of initial state radiation 

For simulation of the initial state radiation (ISR), it is convenient to represent the radia- 
tive correction in the form 



1 + 5' 



71 W 



^ (12) 



,1-/3 ' 



where (3 = (a/rt)(L — 1), r = E~jE h , and E 1 is the energy of the ISR photon. 

The function under the integral can be interpreted as the energy spectrum for photons 
radiated from the initial state. Indeed, at Q 2 = 1 4- 100 GeV 2 the parameter /3 is small 
((3 = 0.033 -j- 0.044), and this function coincides approximately with the energy spectrum 



a 



for hard photons, radiated from the initial state 

dN a(L - 1 



(13) 



dr txt 

For simulation of the extra photon emission, we replace the four- dimensional integration 
in Eq. (JH]) to five-dimensional one with r as the outermost integration variable 



da 



a ( 3 
1 + - -L-l 



7T V 4 



r 



^da B dr (14) 



The vacuum polarization correction is included by the substitution 

a 2 -> a 2 (l + 5 VP (Ql))(l + 5 VP {QD) (15) 

in the Born cross section das- 

In simulation of the initial state radiation, the approximation is used that the photon 
is emitted strictly along the initial direction of the radiating electron. Since the energy 
of the photon is restricted by the condition r < r max , we expect that this approximation 
does not lead to a significant systematics in determination of the detection efficiency. Note 
that selection criteria used in data analysis should provide the fulfillment of the condition 
r < r max for both experimental and simulated events. 

To increase simulation efficiency, the variable r is initially generated according to the 
/3 /r 1-/3 ° distribution with f3 = P(Qmm)> wriere Qmin * s a lower bound on the tagged- 
electron Q 2 for simulated single-tag event. If the generated value of r is higher than a 
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threshold r min , the photon is added to the list of final particles in an event. The scattered 
e + and e~, and the pseudoscalar meson are then generated in the frame with the shifted cm. 
energy of 2Ef,\/l — r. If r < r min , the photon is not generated, and the cm. energy is not 
shifted, but the radiative correction factor in the cross section (see Eq. ( Fl4|) ) is calculated. 



3.2. Simulation of final state radiation 



The final state radiation (FSR) is simulated after the generation of the two-photon event. 
The final electron scattered at a large angle is "decayed" to e + 7 with some probability. The 
final-meson four-momentum is then modified to provide the energy and momentum balance. 
The probability of the emission of the photon with the energy greater than E ljmin equals 



P(Q 2 ,x r 



a 



tt(1 + 5') 



(L - 1) In 



-L 



(16) 



where x m i n = E^ ym i n /E, and E is the electron energy before FSR simulation. This formula 
is obtained by integration of the FSR photon spectrum given by Eq. (23) of Ref. Q. The 
Q 2 dependence of the FSR probability calculated for x min = 0.1, 0.01 and 0.001 is shown in 
Fig. EJ 
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FIG. 5. The Q 2 dependence of the final state radiation probability. 



The photon energy E 1 and angle 9 1 with respect to the electron direction before radiation 
are generated according to the following distribution function J9]: 

diV a 1 - x + x 2 /2 



dxd cos 6~ 



TTX 1 — COS 6~ 



(17) 



where x = E^/E, /3 = \/l — m 2 /E' 2 , and E' is the electron energy after the photon emission. 
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4. COMPARISON WITH OTHER GENERATORS 



The comparison of the total cross sections in the no-tag mode obtained with GGRESRC 



and the two other generators of two-photon events, GGRESPS [6| and TWOGAM ll|, was 
performed. The results of Monte-Carlo calculations are identical for all the three generators, 
if the mass of the meson, its two-photon width, and Q 2 -dependence of the form factor are 
set to be the same in the generators. The GGRESRC and GGRESPS use the same formula 
(Eq. dSJ)), but different orders of integration over the invariant variables. The TWOGAM 
generator was developed for the CLEO measurements of the meson-photon transition form 
factors 12]. It is based on the BGMS formalism [lj (see Eq. d^J)) and uses the completely 
different integration variables, the momenta of the final electrons. 

For GGRESRC in the regime without radiative corrections and TWOGAM, the com- 
parison of the Q 2 spectra, obtained for the process of the 7r° production in the single-tag 
mode, was performed. The spectra was found to be in agreement within the Monte-Carlo 
statistical errors. 



5. GENERATOR PARAMETERS 

The parameters of the event generator are listed in Table IIHI The recommended values 
for the parameters Rmax, Rmin, and Kmin are given in brackets. To simulate no-tag events, 
the parameters QlSmin, Q2Smin, and IRad should be set to zero. The regime with radiative 
correction (IRad=l) is used only in the single-tag mode. 

The resonances decay modes implemented in the generator are listed in Table 4. The 
decay models used are described in Ref. [6j. If parameter IMode equals 0, the meson decay 
is not simulated. 

In general terms the GGRESRC simulation algorithm is the following: 

• The electron and positron collide in the cm. frame (So). In this frame the positive 
z-axis is defined to coincide with the e~ beam direction. 

• The emission of a hard photon from the initial state is simulated. The photon is 
emitted along the collision axis. If the ISR photon energy is greater than r min E b , the 
photon momentum is stored in the list of final particles. 



11 



TABLE III. Parameters of the generator GGRESRC. 



Name 


Description 


Eb 


beam energy (GeV) 


IR 


produced meson: n° (= 1), n (= 2), n' (= 3), ri r {= 4), m, (= 5) 


IMode 


meson decay mode (see Table lIVp 


KVMDM 


form factor model: constant (= 0), VDM (= 1) 


Itag 


tagged particle: positron (= 1), electron (= 2), mix (= 3) 


IRad 


simulation with/without radiative correction (= 1/0) 


Rmax 


maximal energy of ISR photon in units of Eb (0.1) 


Rmin 


minimal energy of ISR photon in units of Eb (10 — 4 ) 


Kmin 


minimal energy of the FSR photon ( 0.001 GeV ) 


QlSmin 


minimal momentum transfer squared to the untagged electron 


QlSmax 


maximal momentum transfer squared to the untagged electron 


Q2Smin 


minimal momentum transfer squared to the tagged electron 


Q2Smax 


maximal momentum transfer squared to the tagged electron 


Fmax 


maximum weight of events 



• The scattered electrons and the resonance are generated in the new cm. frame (Si) 
with the cm. energy 2E b y / l — E^j 'E b ; Si = S if E 1 jE h < r min . 

• The final state radiation is simulated. If the photon energy is greater than k m i n , its 
parameters are stored in the list of final particles. The momenta of the tagged electron 
and the produced meson are modified. 

• The meson decay is simulated. 

• The momenta of the final particles are transformed from Si to So frame. 

When required statistics is collected, the total cross section for the two-photon process with 
radiative corrections is calculated and printed. 
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TABLE IV. The meson decay modes in GGRESRC. If IMode=0, the meson decay is not simulated. 



Meson 


IMode 


Decay channel 






Branching 














fraction [13] (%) 




1 


7T° - 


-4. 2^ 
-t z 7 






98.798 




2 


7T°- 


-> e + e~7 






1.198 




1 


77 — 


► 27 






39.31 


V 


2 


7? — 


► 3ir°, ir° - 


->■ 


2 7 


31.4 




3 


77 — 




7T 


-> 27 


22.457 




4 


77 — 








4.6 




1 


77' — 








2.1 


rf 

1 


2 


'/ 


-7* TT^ 7T~ Tj ^ 


77 


-»■ 27 


17.532 




3 


TV — 


7r^~7r _ 77, 


77 
/ 


-> 3tt°, vr° -)• 27 


14.004 




4 


??' — 


-7* 7T^~7T 7y, 


77 


—J- 7T + 7r~7r , 7T° — > 27 


10.016 




5 




7 


^* /I /I //, 


77 


—7- 7T + 7r~7 


2.0516 




6 


??' - 
'/ 


-> 27T^77. 77 


— > 


2-y vr° — v 2t 


7.943 




7 


'/ 




-> 


3vr°, 7T° -»• 27 


6.3445 




8 


7?' — 


-V 27T^77 77 


-> 


7r + 7T _ 7T , 7T° — > 2^ 


4.5375 




9 


7?' - 


-> 27T°7/, 7/ 


-> 


7r+7r~ 7, 7T° — > 27 


0.9294 




10 


7/ - 


+ A, p° 


-)■ 




29.4 


Vc 


1 


% " 


-> K s K+it- 


+ C.C. 


2.33 




2 




^2 7 






0.024 


Vb 


1 




-> 27 









6. EXAMPLE OF e + e~ -> e+e vr SIMULATION 

In this section, some distributions for the process e + e~ — > e + e~ + 7r°, (tt° — > 27) obtained 
with the generator GGRESRC, are presented. In Table |V] the parameters of the generator 
used in the simulation are listed. 

At these parameter values, 57% of events do not contain extra photons, 22%, 28%, and 
6% of events contain ISR photon, FSR photon, and both ISR and FSR photons, respectively. 
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TABLE V. The simulation parameters used for calculation of the distributions shown in Figs. [61-1131 



Parameter 


Value 


Comment 


Eb 


5.29 


beam energy (GeV) 


IR 


1 


meson: ty® 


IMode 


1 


decay mode: 7r° — > 2"f 


KVMDM 


1 


VDM form factor (Eq. ©) is used 


Itag 


1 


the final positron is tagged 


IRad 


1 


radiative corrections are simulated 


Rmax 


0.1 


maximal energy of ISR photons in units of Eb 


Rmin 


10~ 4 


minimal energy of ISR photons in units of Eb 


Kmin 


0.001 


minimal energy of FSR photons (GeV) 


QlSmin 





Qmin for e" (GeV 2 ) 


QlSmax 


1.5 


Qmax for e~ (GeV 2 ) 


Q2Smin 


1.5 


Q 2 min for e+ (GeV 2 ) 


Q2Smax 


9 


Qliax for e+ (GeV 2 ) 



The obtained cross section of the process and average radiative correction are: a =0.99 pb, 
5 =-0.6%. 

The calculated cross section as a function of the restriction on Q\ (the values of the other 
parameter are equal to those in Table |V]) is shown in Fig. |6j One can see that at Q\ max ~ 
1.5 GeV 2 the cross section reaches an asymptotic value. 

The energy spectra of tagged electrons obtained with and without radiative-correction 
simulation are shown in Fig. [7J It is seen that emission of extra photons significantly changes 
the shape of this spectrum. 

Fig. |8] shows the energy spectra of the photons from the 7r° — > 27 decay in the process 
e + e~ — > e + e~7r°. The energy spectra of photons emitted by the initial and final electrons 
are presented in Fig. [HI 

The polar-angle distributions of tagged electrons is shown in Figs. [10J It is seen that 
at E b = 5.29 GeV the cut Q 2 min > 1.5 GeV 2 corresponds to the minimum scattering 
angle of about 13°. This is in agreement with an estimate for small scattering angles 
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FIG. 6. The e + e — > e + e 7r° cross section as a function of the limit on Q\. 
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FIG. 7. The energy spectra of tagged electrons from the process e + e — > e + e tt° calculated 
without (left panel) and with (right panel) radiative correction simulation. 
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FIG. 8. The energy spectra of photons from the 7r° — > 27 decay in the process e + e 
calculated without (left panel) and with (right panel) radiative correction simulation. 
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9 Q/ (Eb ■ E'Y/2, where energy of the scattered electron E' w E^. 
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FIG. 9. The energy spectrum of ISR (left panel) and FSR (right panel) photons in the simulation 
of the process e + e~ — > e + e~7r°. 
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FIG. 10. The polar-angle distributions of tagged electrons from the process e + e — > e + e ir 
obtained without (left panel) and with (right panel) radiative correction simulation. 

The polar angle distribution of FSR photons is shown in Fig. [TTJ Since the FSR photon is 
emitted predominantly along the tagged-electron direction, the photon angular distribution 
is very close to that for the electron. 

The polar-angle distribution of photons from the tt° — > 2 / y decay is shown in Fig. [12j 
Photons have wide distribution, which becomes more uniform with account of radiative 
corrections. 

In Fig. Unithe distribution of the missing mass in the process e + e~ — > e + e" + tc° is shown. 
The missing mass is calculated as a/G°i + Vi ~~ k — p' 2 ) 2 , i-e. we assume that only the tagged 
electron and the two photon from 7r° decay are detected. The narrow peak at zero mass 
contains events (57% of the total number of events), which do not have extra ISR or FSR 
photons. It is seen that emission of the extra photons leads to significant widening of the 
missing mass distribution. 
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FIG. 11. The polar- angle distributions of FSR photons from the process e + e — > e + e it . 




FIG. 12. The polar-angle distributions of photons from the tt — > 2j decay in the process e + e — > 
e + e~ir° obtained without (left panel) and with (right panel) radiative correction simulation. 

7. PROGRAM COMPONENTS 



7.1. Common blocks 



COMMON /GGRSTA/Sum , Es , Suml , Sum2 , Fm , Fml , NOBR , Nact , Ngt 

Purpose: to collect simulation statistics. 

Sum Real*8 used for calculation of the total cross section 

Es Real*8 used for calculation of the cross-section error 

Suml Real*8 used for calculation of the average form factor 

Sum2 Real*8 used for calculation of the average radiative correction factor 

Fm Real*8 maximum weight of event 

Fml Real*8 maximum weight of event in FSR simulation 

NOBR Integer*4 number of calls of the generator 
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FIG. 13. The missing mass distribution for the simulated e + e —> e + e ir° single-tag events. 



Nact Integer*4 number of the generated events 

Ngt Integer*4 number of events with the weight greater than Fmax (see /GGRPAR/) 

COMMON / GGRPAR/Eb , Rmas , Rwid , Rg , Rm , Fmax , Rmax , Rmin , t 1 imin , t 1 imax , 
t2imin , t2imax , Kmin , Fmaxl , IR, IMode , KVMDM , ITag , IRad 



Purpose: 


simulation parameters. 


Eb 


Real*8 


beam energy (GeV) 


Rmas 


Real*8 


meson mass (GeV) 


Rwid 


Real*8 


meson total width (GeV) 


Rg 


Real*8 


meson two-photon width (keV) 


Rm 


Real*8 


meson mass in the current event (GeV) 


Fmax 


Real*8 


expected maximum weight of event 


Rmax 


Real*8 


maximal energy of ISR photon in Eb units 


Rmin 


Real*8 


minimal energy of ISR photon in Eb units 


tlimin 


Real*8 


minimal value of t\ (GeV 2 ) 


tlimax 


Real*8 


maximal value of t\ (GeV 2 ) 


t2imin 


Real*8 


minimal value of t% (GeV 2 ) 


t 2 imax 


Real*8 


maximal value of t2 (GeV 2 ) 


Kmin 


Real*8 


minimal energy of the FSR photon (GeV) 


Fmaxl 


Real*8 


maximum expected weight for FSR simulation 


IR 


Integer*4 


meson type 



Imode Integer*4 meson decay mode 

KVMDM Integer*4 form factor model 

ITag Integer*4 tagged particle 

IRad Integer*4 switch for radiative correction calculation 



COMMON /GGRCON/Alpha , PI , EM , mPiO , mPi , mEta , mEt ap , mKs , mKc , mRho , m Jps 
mUps,BrPiO(2) ,BrEta(4) ,BrEtaPrim(4) ,BrRho,BrTot 



Purpose: constants. 




Alpha 


Real*8 


fine structure constant (1/137.03604) 


Pi 


Real*8 


7T (3.14159265) 


Em 


Real*8 


electron mass (0.00051099891 GeV) 


mPiO 


Real*8 


7T° mass (0.1349766 GeV) 


mPi 


Real*8 


tt ± mass (0.13957018 GeV) 


mEta 


Real*8 


r] mass (0.547853 GeV) 


mEtap 


Real*8 


i mass (0.95766 GeV) 


mKs 


Real*8 


K s mass (0.497614 GeV) 


mKc 


Real*8 


X ± mass (0.493677 GeV) 


mRho 


Real*8 


p° mass (0.77549 GeV) 


mJpsi 


Real*8 


J/i) mass (3.096916 GeV) 


mUps 


Real*8 


T mass (9.4603 GeV) 


BrPiO(2) 


Real*8 


7T° decay branching fractions 


BrEta(4) 


Real*8 


rj decay branching fractions 


BrEtaPrim(4) 


Real*8 


r]' decay branching fractions 


BrRho 


Real*8 


branching fraction of the decay p° — > tt + tt 


BrTot 


Real*8 


total probability of the decay chain 



COMMON /GGREV/pPart(4,25) ,mPart(25) , Type (25) , Mother (25) ,Npart 
Purpose: final particle parameters (up to 25 particles). 
pPart(l-3,i) Real*8 momentum of i-th particle (GeV) 
pPart(4,i) Real*8 energy of i-th particle (GeV) 
mPart(i) Integer*4 mass of i-th particle (GeV) 

Typed) Integer*4 type of i-th particle 
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Mother (i) Integer*4 index of parent of i-th particle in /GGREV/ 
Npart Integer*4 total number of particles in /GGREV/ 

In the common block /GGREV/: 1-st and 2-nd particles are the scattered electrons, 3-rd 
particle is the produced resonance, 4-th e.t.c. particles are the ISR photon (if exists), the 
FSR photon (if exists), resonance decay products. 

COMMON /GGRPOL/SETS (7330) , SETPOL (7330) 
Purpose: vacuum polarization correction. 

SETS Real*8 momentum transfer squared (GeV 2 ) 

SETPOL Real*8 value of the vacuum polarization correction 

Common blocks for internal use: /GGRARIP/, /GGRFUC/. 



7.2. Subroutines of the generator 



GGRESRC 


the main subroutine 


GGRDEC2G 


simulation of resonance decay to 27 


GGRESEND 


print out of simulation results 


GGRESINI 


initialization 


GGRETCD 


simulation of r\ c decays 


GGRETD 


simulation of rj decays 


GGRET1D 


simulation of 7/ decays 


GGRET1D1 


simulation of the decays rj' — > 7i + 7i^i] and ii 7r r) 


GGRET1D2 


simulation of the decay rf — > p°7 


GGRFSR 


FSR simulation 


GGRFVP 


filling the common block /GGRPOL/ 


GGRINV 


simulation of the invariants t 2 , t±, s 1 , s 2 


GGRLMOM 


calculation of the laboratory momenta of the final electrons and meson 


GGRLOR 


Lorentz transformation 


GGRPIOD 


simulation of ir° decays 


GGRPI0D1 


simulation of the 7r° — > e + e~7 decay 


GGRPREV 


print out of one event 


GGRRNDM 


wrapper of a pseudo-random numbers generator 


GGRSPC3 


simulation of the three particle phase space 
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7.3. Double-precision functions 



GGRFU 



GGRPOLAR 



calculation of the vacuum polarization correction 
function used by the subroutine GGRFSR 



GGRFVMDM 



calculation of the form factor in the vector dominance model 



7.4. Library subroutines 



In the generator we use following functions from the CERN program library: 
RANLUX generation of pseudo-random numbers uniformly distributed in the interval (0,1); 
DZEROX computing a zero of a real- valued function f(x) in the given interval [a, b]. 



The event generator GGRESRC for simulation of the two-photon process e + e~ — > e + e~R } 
where R is a pseudoscalar meson, has been developed. The generator allows to efficiently 
generate two-photon events in the single-tag mode, when one of the final electrons is scattered 
at a large angle and may be detected. In this mode simulation of radiative corrections has 
been implemented in the generator including extra photon emission from the initial and final 
states. 

The generator is used for simulation of experiments with the BABAR detector on mea- 
surements of the photon- meson transition form factors (see, for example, Refs. ) , and 
for simulation of two-photon experiments with the KEDR detector at VEPP-4M collider. 



The work is partially supported by the RF Presidential Grant for 
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